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Abstract. We pres ent an extinction m ap of the Polaris molecular cirrus cloud derived from star counts and 



compare it with the Schlegel et al. (1990 extinction map derived from the far-infrared dust opacity. We find that. 



within the Polaris cloud, the Schlegel et al. (1998) Av values are a factor 2 to 3 higher than the star count values. 
We propose that this discrepancy results from a difference in rpm/Av between the diffuse atomic medium and 
the Polaris cloud. We use the difference in spectral energy distribution, warm for the diffuse atomic medium, cold 
for the Polaris cloud, to separate their respective contribution to the line of sight inte grated infrared emissi on and 
find that the tfir/Av of cold dust in Polaris is on average 4 times higher than the [Schlegel et al. (1998| ) value 



for dust in atomic cirrus. This change in dust property could be interpreted by a growth of fluffy particles within 
low opacity molecular cirrus clouds such as Polaris. Our work suggests that variations in dust emissivity must be 
taken into account to estimate Av from dust emission wherever cold infrared emission is present (i.e. molecular 
clouds). 

Key words. ISM : clouds - ISM : dust, extinction - ISM : individual object: Polaris Flare - Infrared: ISM: 
continuum 



1. Introduction 

The infrared sky images provided by the Infrared 
Astronomy Satellite (IRAS) and the Cosmic Background 
' Explorer (COBE) have considerably improved our knowl- 
' edge of interstellar dust and of the spatial distribution of 
interstellar gas at high Galactic latitude. With the IRAS 
data it became clear that interstellar dust comprises small 
particles stochastically heated by the absorption of pho- 
tons to temperatures higher than the equilibrium temper- 



ature of large grains (Desert et al., 1990). Far from heating 



sources, these small particles make the 12, 25 and a sig- 
nificant fraction of the 60 /im IRAS emission. Only the 
100 fim emission comes from large grains emitting at the 
equilib rium temperature sets by the balance between heat- 



particular, [Schlegel et al. (1998|) (hereafter SFD98) com- 
bined the angular resolution of IRAS 100 /xm data and 
the wavelength coverage of DIRBE to make a map of the 
dust far-infrared (FIR) opacity. They have shown that, 
in regions of low extinction {Ay < 0.3) at high Galactic 
latitude, the FIR dust optical depth, tfir, is tightly cor- 
related with visible extinction. By extrapolating this cor- 
relation to the whole sky they produced an all-sky map of 
visible extinction which is being used for a wide range of 
purposes. 

The IRAS and COBE observations show that the spec- 
tral energy distribution of dust emission varies in the 
ISM and in particular from the diffuse atomic medium 



ing an( |l cooling. COBE extended the IRAS observations 



to molecular clouds (e.g. Laureijs et al., 1991; Aberge] 
et al., 1994; Lagache et al., 1998| ). These variations have 



to the far-infrared/sub-millimeter emission. These data 
have been used to measure large grain properties (tern- 



terstellar matter ( 


Wright et al., 1991; 


Reach et al., 1995; 


Boulanger et al., 1996 


; Lagache et al., 1998; 


Arendt et al., 


199S; schlegel et al., 1998 




Finkbeiner et al., 199E). In 
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been interpreted in terms of changes in dust composition, 
in particular variations in the abundance of small grains. 
What is still a matter of debate is whether they also im- 
ply changes in the properties of large grains. In this paper 
we address this question by looking for variations in the 
tfir/Av ratio. In a first part of the paper, we compare 
the SFD98 extinction map with an independent estimate 
of Ay derived from star counts over the Polaris Flare, a 
high latitude cirrus with CO emission and an extinction 
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around 1 mae dHeithauscn fc Thaddeus, 1990| ; PBernard 
et al., 1^9991 ). We find that the SFD98 extinction IS on av- 
erage twice that estimated from star counts. In the second 
part, we investigate various explanations of the extinc- 
tion discrepancy and propose to relate it to a difference 



in TFrn/^v between dust in low extinction cirrus and the 



colder dust associated with the Polaris Flare. 



2. Extinction in the Polaris Flare 

2.1. Extinction map from star counts 



We have used the USNO-PMM ( [Monet, 1996|) B pho- 
tometry and the method described in Cambresy (1999) 
to build a visual extinction map (assuming Ry = 3.1) of 
the Polaris cirrus cloud (Fig. |l|). The method consists of 
star counts with variable resolution in which the local stel- 
lar density is estimated with the 20 nearest stars of each 
position. The extinction is derived from the star density 
(D) as follows: 



Ab 
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where a is the slope of the B luminosity function and Drof 
the reference stellar density in the absence of extinction. 
The high galactic latitude of the Polaris Flare [h ~ 27°) 
prevents the contamination by other clouds on the same 
line of sight. The upper lim it to its distance of about 240 
pc ( Heithauscn et al., 1995 ) ensures that the cloud is close 
enough to derive the extinction from star counts without 
significant contamination by foreground stars. 

The extinction map (Fig. |^) is obtained by taking for 
log-Di.cf(&) the value derived from the linear fit of log I? 
versus the galactic latitude h outside the Polaris Flare, 
i.e. logi^rcf = P + a X b. This corrects the stellar den- 
sity variations due to variations in the length of the line 
of sight through the stellar disk, but ignores extinction 
from diffuse interstellar matter outside the Polaris cloud. 
To account for this diffuse uniform or slowly varying ex- 
tinction which is not seen by star counts, we measured 
the SFD98 map variations with respect to the galactic 
latitude for regions with {Ay < 0.2). For b ranging from 
25° to 34° and I ranging from 105° to 143° we find that 
Ay(SFD98) - Ay(starcount) w 0.20 - 7.3 lO'^ x b. We 
have used this relation to correct the star count extinction 
values in order to make the comparison with the SFD98 
map; the absolute value of the resulting correction is lower 
than 0.05 mag. After this correction, both our extinction 
map and the SFD98 map measure the total extinction. 
This is necessary to compare each map together and with 
DIRBE far-infrared data. Extinction values outside the 
Polaris cloud (where Ay < 0.2) match those of the SFD98 
map. For 17 < 6 < 25 the linear relation we measured is 
no longer valid and it becomes hard to make such a cor- 
rection because of the lack of areas with low extinction in 
the SFD98 map. Since we are interested only in directions 



represented by white crosses which are all at & > 25 in 
Fig. |l], there is no need to try to include diffuse extinction 
for 6 < 25. The star count Polaris extinction map agrees 
with extinction values derived from color excesses based 
on CCD images of a 1 deg^ section of Polaris (Zagury 
et al., 1999| ). 

The resolution of the resulting map is about 8'. 
Statistical uncertainties in the determination of the ex- 
tinction come from star counts. The distribution of stars 
follows a Poisson law with a sigma equal to ^/N, where 
N is the number of counted stars. For iV = 20 we obtain 
AAy = ±o'23- Moreover, the Ry parameter introduces 
a systematic uncertainty. In the absence of any measure- 
ment, we have assumed Ry = 3.1 which corresponds to 
the average value in the diffuse interstellar medium. A 
higher value of 5.5, that is common but not systematic 
in dense cores, would multiply our extinction values by a 
factor 1.12. 
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Fig. 1. Visual extinction map of the Polaris Flare from 
B star counts using the USNO-PMM catalogue. White 
crosses represent the DIRBE pixels used to compare the 
star count based extinction map with the SFD98 extinc- 
tion map based on FIR emission 



2.2. Comparison with SFD98 

In Fig. H, the ratio between SFD98 and star count ex- 
tinction values is plotted against the dust temperature as 
derived by SFD98. For this we have smoothed both ex- 
tinction maps to the DIRBE resolution. The DIRBE pix- 
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for this comparison are marked with white crosses fc Guhathakurta (1999| ) have presented a study of cir- 



els usei 

in Fig. 111. The SFD98 Ay values are systematically higher 
than those derived from star counts. The extinction ratio 
decreases for increasing dust temperatures. The mean ra- 
tio between the SFD98 and star count extinction values is 
2.1. Fig. ^ shows the spatial distribution of the difference 
between SFD98 and our extinction map. 



rus clouds with U BVRI photometries and have com- 
pared their extinction derived from star counts, color ex- 
cess, UBV color-color diagrams and BVI color-color di- 
agrams. They found for their highest extinction cirrus, 
in Corona Australis, a discrepancy with the SFD98 map 



comparable to that found here for the Polaris Flare. Arcs 



& Goodman (1999) have compared the extinction toward 
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Fig. 2. Comparison of the SFD98 extinction with star 
count extinction in the Polaris Flare versus the SFD98 
temperature 
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Fig. 3. Difference between SFD98 and the star count ex- 
tinction map at the DIRBE resolution. The two maps are 
in agreement outside the molecular cloud 

The Polaris star counts raise a question about the va- 
lidity of the SFD98 extinction map in regions of higher 
extinction that in those where they calibrated their con- 
version from FIR dust opacity to Ay. Previou s studies 
have already pointed out similar discrepancies. Szomoru 



the Taurus cloud using four different techniques and they 
also concluded that the SFD98 map overestimates the ex- 
tin ction by a factor of 1.3-1.5 . These results are confirmed 
by von Braun fc Mateo (2001 ) who compared a color study 
of the globular cluster NGC 3201 (6 = 8.6°) with the 
SFD98 extinction map for this line of sight. They con- 
cluded their analysis on a reddening overestimation in the 
SFD98 map. 



2.3. Discussion 

The extinction data presented in this paper together with 
similar data on other molecular clouds all lead to the same 
conclusion: the SFD98 extinctions are larger than the val- 
ues derived from star counts in molecular clouds where 
Av>l . 

Cloud structure on angular scales smaller than the res- 
olution of the star counts can make the measured Ay 
smaller than the mean value of Ay - The relation between 
the mean extinction and the stellar density is not linear 
but logarithmic (Eq. |l|) whereas the FIR integration over 
a cell is linear. This could thus explain the extinction dis- 
crepancy. Rossano (1980| ) has quantified this effect. Using 
his results, we found that a difference of a factor 2 for a 



visual extinction of 2 magnitudes would require a surface 
filling factor lower than 0.1. Such a low surface filling fac- 
tor is incompatible with studies of the dumpiness of dust 
extinction which conclude on a smooth distribution of the 



19991 ) 



dust in molecular clouds ( Thoravalet al., 1997 ; Lada ct al 



Another plausible explanation would be a shortcoming 
in the SFD98 dust temperature determination in molec- 
ular clouds which ignores temperature variation. This 
should be considered due to the coarse resolution of the 
SFD98 temperature map (~ 1°) compared to the dumpi- 
ness of the brightest structures in cirrus clouds. For exam- 
ple, for a given FIR brightness, an underestimation of the 
dust temperature translates into an overestimation of the 
dust opacity. Observations of molecular clouds at higher 
angular resolution than DIRBE do show temperatures lo- 
cally lower than those of SFD98. For example, PRONAOS 
observations of a piece of the Polaris cirrus gave a dust 



temperature of 13.5 K (Bernard et al., 1999), significantly 



lower than the value of 16 K at the same position in the 
SFD98 temperature map. This example clearly shows that 
the SFD98 temperature is an effective temperature which 
represents a mean intensity value weighted over the dust 
seen within the solid angle of their study. The impact of 
the beam (and line-of-sight) averaging on the SFD98 tem- 
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peratures and opacities can be discussed on the basis of 
a simple model where we assume that the emission comes 
from a mixture of two distinct components. In this cal- 
culation, we assume an emissivity law (oc ly^) and fit a 
single black body to a combination of two black bodies 
at different temperatures. The nffnrtive optical depth is 



3.1. Evidence for distinct emission components 

The IRAS data have shown that no emission from aro- 
matic hydrocarbons and very small grains at 12, 25 and 
60 fim is seen towards the dense gas traced by the mil- 
limeter transitions of ^^CO (Laureijs et al., 1991; Aberge] 



aJwaysl smaller tha.n the sum of the optica,] depth for the 
two temperatures. Thus if the tfir/Av is the same for 
both dust components, our simple model shows that the 
SFD98 method always lead to an underestimation of Ay, 
the reverse of the discrepancy observed with star counts 
(Fig. More generally, when a range of temperatures 
is present within the beam, the effective FIR dust opac- 
ity that is measured by SFD98 is always smaller than the 
mean opacity. Temperature variations within the beam 
thus cannot explain the extinction discrepancy between 
SFD98 and star counts. 

The SFD98 map is a tfir map and not an Ay map. To 
go from one to the other one needs to know the tfir/Av 
ratio. We propose to explain the extinction discrepancy by 
variations in the tfir/^v' ratio between the low extinction 
regions used for the calibration of the Ay and the Polaris 
molecular cirrus. Variations in the Ry ratio from diffuse 
to dense clouds are interpreted as evidence for variations 
in the o ptical properties of la rge dust grains in molecular 
clouds ( ICardelli et al., 1989 ). In our interpretation, the 
extinction discrepancy would thus be an additional signa- 
ture ofi thc evolution of dust from diffuse to dense clouds. 



et al., 1994). In practice, the molecular emission from this 



dense gas is observed to be well correlated with the dif- 
ference between the 100 and 60 /xm IRAS brightnesses. 



Lagache et al. (1998D combined the FIR DIRBE bands 
to the 60 /im data to separate along each line of sight 
the emission associated with the diffuse ISM and dense 
gas where both are present. With the long wavelength 
bands of DIRBE (100, 140 and 240 /xm), they were able 
to determine a dust temperature for each of these two 
emission components. Their results show that the dust 
associated with the dense g traced by the difference 
between the 100 and 60 /im brightness, is systematically 
colder (T ~ 15 K) than dust in the diffuse atomic in- 
terstellar medium (T ~ 17.5 K). We refer to these two 
emission components as cold and warm. The distinction 
between cold and warm emission components character- 
ized by different dust temperatures and the abundance of 
small grains emission is corroborated by sub-mm obser- 
vations at higher angular resolution for which the spatial 
separation of the infrared emission from diffuse and dense 
gas is easier to distinguish than with the DIRBE data 
(Bernard et al., 1999; Stepnik et al., 2001; Laureijs et al 



1996) 
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We propose to relate the extinction discrepancy to 
a variation in the Tpm/Ay ratio between the cold and 
warm dust components. In their analysis of the DIRBE 
data, SFD98 did not consider the possible presence of dis- 
tinct emission components with different dust tempera- 
ture along the line of sight. They assume that large grain 
properties are everywhere the same and consequently that 
the temperature variati ons are exclusively du e to changes 
in dust heating. In the Lagache et al. (1998 ) description 
of the data, the observed variations in the dust effective 
emission temperature are to a large extent due to vari- 
ous degrees of mixing between the warm and cold com- 
ponents. In a second paper, Finkbeiner et al. (199£ ), the 
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Fig. 4. Effective to total optical depth ratio versus warm 
component fraction. Cold and warm components are de- 
scribed by the modified Planck function P = B^{T) x 
(A/Ao)-2 with Tc = 13.5 K and = 17.5 K, respectively. 
Optical depths are derived by fitting each spectrum with 
a single modified Planck function. 



3. Warm and cold dust components 



same authors as SFD98, showed the necessity of includ- 
ing a very cold emission component with T 9K in their 
analysis to improve the fitting of the FIRAS data at long 
wavelengths. However, this very cold component and the 
cold component of Lagache et al. (1998) are distinct (Note 
that they are not mutually exclusive and that both may 
be required to correctly describe the data) . The cold com- 
ponent, derived from DIRBE data, is spatially correlated 
with known molecular clouds; the very cold component, 
deduced from FIRAS long-wavelength spectra is, on the 
contrary, present in all directions. For the subject of this 
paper the key difference betw een both components is that 
in the Finkbeiner et al. (1999| ) analysis the very cold com- 
ponent accounts for a fixed fraction of the emission, inde- 
pendent of sky position (fixed opacity ratios between the 
very cold and warm components, with temperature given 
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by a power law) while in the Lagache et al. (1998) anal- 
ysis the cold fraction varies spatially and is zero in the 
high latitude regions used by SFD98 for the calibration 
of the Tpiji/Av ratio. Due to this spatial dependence, a 
difference in dust properties between the cold and warm 
emission components, namely in the Tpm/Ay ratio, will 
limit the domain of the SFD98 conversion of tfir into 
Ay and could account for the observed discrepancy with 
star counts. This statement is quantified in the following 
sections. 

3.2. Separation of the warm and cold components 

A separation of the warm and cold contributions to 
the infrared brightness is necessary in order to derive a 
temperature and FIR opacity for each component. To 
do this, we apply a method similar to that of Lagache 
et al. (199^ ) to DIRBE images of the Polaris Flare. 
Let /(A) (A=100, 140, 240 ^m) be the DIRBE emission 
at 100, 140 and 240 /im respectively, and i?(A, 60) the 
/i/(A)//i/(60) flux ratio observed in cirrus clouds. The cold 
maps are computed at each wavelength according to the 
relationship: /(A)c = /(A) - i?(A, 60) x /(60). 

We first determine i?(100, 60) using the correlation 
diagram of the diffuse cirrus emission around the Polaris 
cloud. We obtain i?(100, 60) = 4.0 ± 0.3. At 140 and 
240 /im the correlation diagrams are much more noisy. 
Therefore we prefer to compute i?(140, 60) and i?(240, 60) 
using i?(100, 60) and assuming a large grain temperature 
of 17.5±0.5 K. We thus obtain: i?(140,60) = 7.9 ± 1.1 
and i?(240,60) = 6.4 ± 1.4. 

The produced cold maps have non-zero emission at a 
large scale due to 1) zodiacal residual emission (mainly 
coming from the 60 /im map) and 2) non-correction of 
the Cosmic Infrared Background. Therefore, under the 
assumption that the cold dust is distributed with limited 
angular extent in molecular clouds, we remove the low fre- 
quency structures using a 18° x 18° median filter. We thus 
obtain, at the end, maps of cold dust emission at 100, 140, 
and 240 /im which are clumped on scales smaller than our 
filter size. 

Map of statistical uncertainties for the cold emission 
are computed using the DIRBE release error maps: 



the statistical and systematic 100, 140 and 240 /im error 
maps of cold emission, respectively. Temperatures and 
optical depths are derived using the fitting and assum- 
ing a FIR dust emissivity index of 2. Uncertainties given 
by the fitting correspond to the 68.3% confidence level. 

To determine the warm optical depth, we adopt the 
following method. We first remove to the 100 /im map 
the cold 10 /^m emission and th e cosmic FIR background 
value from Lagache et al. (2000 ) . We then assume that the 
warm component has a temperature of 17.5±0.5 K with 
a FIR dust emissivity index of 2. Optical depth is thus 
directly proportional to the 100 /im warm emission. 

The error on the warm emission, including both the 
systematic and statistical errors, is: 



A/, (100), 



A/,(100)2 + A/,(100)2 + A/, (100)2 jB (5) 



Optical depth uncertainties correspond to the minimum 
and maximum optical depth value allowed by the combi- 
nation of A/^(100)w and the 0.5 K error on the assumed 
dust temperature. 

3.3. From far-infrared optical depth to visible 
extinction 

Within the two components model, the total extinction 
can be written as the sum of the extinction in the warm 
and in the cold components: 



Av = A 



V 



Al 



and we have: 



IW _ w 

— ^100 



Al 



( Il2£^ ^ 



\Ax 



T 100 



(6) 



(7) 



(8) 



where the ratio (tioo/^v')w,c are the emissivities of the 
dust at 100 /im for the warm and cold components, re- 
spectively. We take for the FIR-to-visible opacity ratio of 
the warm component, (Tioo/^y)w, the ratio measured by 
SFD98 in low extinction regions. Then, we can derive Ay 
using Eq. ^ and the emissivity of the cold component is: 

( IMHI^ ^ '^loo 
V Av i Av - A^ 



(9) 



A,ta/.(A)c = v/A/.(A)2 + i?(A, 60)2 X AI,{60y 
Systematic uncertainties are derived following: 
Asys/.(A)c = Ai?(A,60) X /,(60) 



(3) 



(4) 



Temperatures and optical depths are derived only for pix- 
els containing significant cold emission (see Fig. |^), i.e. 
cold intensities at 100, 140 and 240 /im greater than 3cr, 
a being estimated at each wavelength using the width of 
the histograms of the cold emission maps. 

Statistical and systematic uncertainties on the cold 
dust temperature and optical depth are computed using 



3.4. Far-infrared emissivity of the warm component 

The dust emissivity for the warm component is taken 
from SFD98. In their paper, they provide the value of 
p ~ Eb-v/Iioo — 0.0184, where /loo is the 100 /im bright- 
ness expressed in MJy sr~^ corrected from the zodiacal 
emission and scaled with the measured dust opacity from 
the pixel dependent dust temperature to a fixed tempera- 
ture of 18.2 K. With this scaling, the conversion from /loo 
to FIR opacity becomes pixel independent. We have: 

Tioo _ -^100 _ 1 

Av ~ Ayx Sioo(r) " Rvp Bioo{l8.2K) 



(10) 
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For the mean Solar Neighborhood extinction law we have 
Rv = Av/Eb-v — 3.1. Moreover, a color correction 
should be applied to the optical depth because FIR fluxes 
are expressed for the effective filter wavelength assuming 
a law. The correction consists of dividing the infrared 
flux by a color factor /C defined in the DIRBE explanatory 
supplement as follow: 



V-""" /real 



(11) 



dv 



where Wn is the filter profile. The value given in the ex- 
planatory supplement is /Cioo(18.2iir) = 0.91125. 

Finally, the emissivity for the warm component is: 



/noo\ /Cioo(18.2K) _^3^^Q-7 

\Av)^ Rvp Bioo{18.2K) 

3.5. Far-infrared emissivity of the cold component 
3.5.1. Result 



(12) 



Following the procedure described in Sect. 3.2 and 3.3 we 
use the extinction map (Fig. [^) and the decomposition 
of the FIR flux in two components. We use the DIRBE 
resolution (the extinction map has been convolved by the 
DIRBE beam) and only pixels with high signal-to-noise 
ratio for both components. Fig. ^ show the temperature 
and the optical depth normalized to lcm~^ for the two 
components. The temperature in the Polaris Flare varies 
from 13 K to 15.5 K. 
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Fig. 6. Emissivity variation versus temperature in the 
Polaris Flare, with the typical statistical (cross in the up- 
per right corner) and systematic (squares) uncertainties. 
The filled square represents the SFD98 value which corre- 
sponds to the warm component emissivity. Dust emissivity 
of the cold component is 4.0 times larger than the warm 
one 



Using Eq. ^ we obtain the values of the emissivity of 

the coW component, ^^^^ i presented in the Fig. |[ These 

values are all above the SFD98 value for the warm com- 
ponent (filled square in Fig. |^). The median value for the 
cold component emissivity at 100 /im is: 

l^lrn^ =5.310-' (13) 

This is 4.0 times larger than the SFD98 value for the 
warm component. 

3.5.2. Uncertainties 

In this sub-section, we compare uncertainties on the dust 
parameters with the dispersion of the points in Fig. ^. For 
the warm component, we are not able to separate the sys- 
tematic from the statistical errors (see Eq. ||) . Since these 
errors are small compared to those on the cold compo- 
nent parameters (by a factor of ~ 3) we focus on this last 
contribution to the observed scatter in Fig. 

Systematic errors. In Sect. we describe errors associ- 
ated with the separation of the FIR flux in two compo- 
nents. Systematic errors on tiqo/Av are dominated by un- 
certainties on i?(A, 60) which are used to derive the emis- 
sion of the cold dust. These systematic errors would shift 
the whole set of data points. In Fig. ^ we give the direc- 
tion and amplitude of these shifts for the median value. 
Extreme values for (^^^ arc 3.610"'^ and 7.910""^ and 

correspond to 2.7 and 6.0 respectively for the ratio of the 
emissivity between the two components. For Ry = 5.5, 
we have already mentioned that star count extinction is 
multiplied by 1.12; the resulting median emissivity ratio 
would then be 3.1 (rather than 4.0). 

Statistical errors. The dispersion of the values represented 
in Fig. ^ suggests an increase of the FIR emissivity of 
the dust when the temperature decreases. However we ex- 
amine the possibility of this trend being linked with sta- 
tistical errors. In Fig. ^ the statistical errors are repre- 
sented by error bars in the upper right corner. Since the 
optical depth is related to the FIR emission by the rela- 
tion Tioo = -^loo/'Sioo (r) , errors on the temperature and 
on the optical depth are correlated: temperature overesti- 
mation implies optical depth underestimation. This effect 
contributes to the observed dispersion in Fig. ^ but can- 
not fully explain its amplitude. We believe therefore that 
an intrinsic dispersion of the tiqq/Av exists and that it 
may be related to an evolution of the emissivity with the 
temperature. 

4. Discussion 

4.1. Dust evolution from warm to cold component 



The submilli meter analysis of the Polaris Flare by Bernarcj 
|et al. (1999 ) has revealed an unexpected low tempera- 
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Fig. 5. Temperature and optical depth normalized at 1cm ^ versus the galactic latitude for warm and cold components. 
Error bars represent the mean total uncertainty (statistical plus systematic) for each component 



ture in this cirrus. They have used the extinction map 
of Fig. |l| to model the effect of the attenuation of the 
interstellar ultraviolet and visible radiation field on the 
dust temperature. They concluded that this cannot ex- 
plain the observed temperature difference between Polaris 
and warm cirrus clouds. They thus suggest that the tem- 
perature difference is related to a change in optical dust 
properties. For a larger FIR to visible/UV extinction ratio, 
grains are able to radiate more efficiently the energy ab- 
sorbed and thus their equilibrium temperature is reduced. 
Thereby, our interpretation of the extinction data is in 
agreement with Bernard et al. (1999| ) conclusions. A simi- 
lar conclusion about the FIR dust emissivity is reached by 
the modeling of PRONAOS FlR/sub-mm observations of 
a molecular filament in Taurus ( ^tepnik et al., 2001 ). In 
this modeling, a change in the tfir/^v ratio between the 
cold dust within the filament and the warm dust outside 
by a factor 3.4 is required to reproduce the brightness pro- 
files across the filament at the various wavelengths. This 
study, unlike, ours does not rely on the empirical separa- 
tion of warm and cold components. 



density. Within the Rayleigh limit (27ra ^ A), Qcxt/a 
is independent of the grain size. For particles large rel- 
ative to the wavelength (27ra ^ A), it is Qoxt which is 
roughly size- independent. In practice, the Rayleigh limit 
applies to all grains in the FIR while the second limit- 
ing case applies to large interstellar grains (a ^ O.l/xm) 
in the UV. Within these limiting cases, one finds that 
TpiR scales as 1/p while tuv scales as 1/ap. For increas- 
ing grain porosity, the effective grain density p decreases 
and tfir increases. Thus, the ratio tfir/tuv scales as the 
grain size a. Qualitatively this also applies to the ratio 
Tpm/Av- (2) The second, more subtle, effect is that the 
optical properties of the composite grains differ from the 
optical properties of their constituents: grain properties 
depend not only on their size but also on their composi- 
tion and structure. In particular, the wavelength depen- 
dence of Qcxt I o. changes with grain composition. For the 
specific case of composite carbon-silicate grains, this leads 
to a significant enhancemen t of the opac ity in the FIR rel- 
ative to that in the visible ( Dwek, 1997 ). Dust properties 
of porous com posite grains have also been quantified nu- 



Th( |; most straightforward explanation for the change 1989|) . In their work, [Mathis fc Whiffen (1989| ) show that 



in dust FIR emissivity between diffuse and molecular 
clouds is grain growth through grain-grain coagulation 
and accretion of gas species (e.g. Draine, 1985). These 
two processes should lead to composite grains with sig- 
nificant porosity. The effect of the formation of compos- 
ite fluffy grains on the dust opacity is well explained by 
Dwek (1997| ). Two effects contribute to an increase in the 
dust opacity per unit dust mass. (1) The first contribu- 
tion results from an increase in the effective grain sizes. 
For spherical grains of size a: Text = T^a^QcxtNd where 
Nd is the dust column density and Qoxt the extinction 
efficiency, and Aid ~ 4/37ra''piVd, where p is the grain 



merically (e.g. BazeU fc Dwek, 1990; Mathis fc Whiffen, 



it is possible to fit the extinction curve and its variations 
with Ry with a size distribution of composite porous car- 
bon-|-silicate grains. Their calculations show that enhance- 
ments of the TpiR/yly ratio by a factor of at least 3 can 
be obtained within the constraints set by the UV to near- 
infrared extinction curve. 

The proposed interpretation of the variations in FIR 
dust emissivity implies that grain coagulation would be 
effective even within cirrus clouds and not only in dark 
clouds and proto-stellar condensations. Observed varia- 
tions in the Ry ratio have also been interpreted as evi- 
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dence for grain growth from diffuse to dense clouds (e.g. 
Kim & Martin, 1996). It will thus be interesting to look 
for a plausible correlation between the changes in rpiR /Ay 
and variations in the visible extinction curve. Note how- 
ever that these two signatures of dust evolution from the 



diffuse ISM to molecular clouds arc not necessarily corre- 



lated because Ry and tfir/Av do not share the same 
dependence on the dust porosity and size distribution 



( pVIathis fc Whiffen, 1989| ). 

In the absence of multiband visible photometric data, 
we have no estimate of Ry in the Polaris cloud and one 



can only make rather qualitative comments. Cardelli et al 
(1989|2have shown that the variations in the extinction 
curve from the UV to the near-UV, to a good approxi- 
mation, can be related to one single parameter, Ry. The 
i?V'-dependent extinction curves all converge in the near- 
IR but the change in the Tpm/Ay ratio proposed to ex- 
plain the extinction data is larger that one would expect 
from the simple extrapolation of this convergence to the 
FIR. For the Cardelli et al. relation, the rpiR/Ay varies 
by a factor 1.2 for a large change in Ry from 3.1 to 5.5. 
If we add the effect of such an Ry change on the star 
count Ay estimates (which needs to be multiplied by a 
factor 1.12 due to the fact that our star counts are based 
on a i? image, see section 2.1), the extinction discrepancy 
between SFD98 and star counts, a factor of about 2, is 
reduced by only a factor 1.34. This paper thus suggests 
that the i?v'-dependent extinction curves should separate 
again in the FIR. 

4.2. Validity of the SFD98 extinction map 

We chose the Polaris Flare for its brightness in the cold 
emission map presented in Fig. |^. The white regions in 
this figure are the regions with cold emission, outside the 
Galactic plane, where the present work questions the va- 
hdity of the SFD98. 



Polaris 



^Taurus 



\ IK ^ i Jgr 



Fig. 7. Cold component o f the dust emission at 240 jim 
from Lagache et al. (1998 ) 



As mentioned in Sect. ^.2l additional studies simi- 



lar to the present work support our results ( 


Szomoru & 


Guhathakurta, 1999; 


Arce & Goodman, 1999; 


von Braun 



fc Mateo, 2001 ) for pieces of the Corona Australis and 



sion in Fig. 0. Since the Szomoru fc Guhathakurta (1999) 
work is a prelude to what can be done systematically 
with the Sloan Digitized Sky Survey (SDSS), it is crucial 
for further works to be aware that the dust emissivity 



varies in the FIR, even for cirrus cloud. Arce & Goodman 
(1999D suggest that the discrepancy results from a call 



Taurus clouds. Both regions have an important cold emis- 



bration bias. In the present paper we show that there is 
no general calibration that could be used in all directions. 
Optical properties of grains vary and the conversion from 
FIR fluxes to opacity cannot be limited to a single scaling 
factor. 

For low galactic latitude (|6| < 10°) it is obvious that 
the SFD98 map has to be used with caution since very 
different physical environm ents are mixed along the line 
of sight. Chen et al. (1999 ) propose a 3-dimensional ex- 
tinction model to constrain the galactic structure based 
on the SFD98 map used toward globular clusters as close 
as 3° to the Galactic center. They found a recalibration 
was better to fit their data and they used the new cali- 
bration for the whole sky. The Galactic plane and even 
more, the Galactic center, are very complicated regions: 
the presence along the line of sight of massive stars which 
heat the interstellar dust and molecular clouds with cold 
dust makes impossible the separation of the different com- 
ponents which contribute to the FIR flux. It is not clear 
how their empirical calibration of the FIR emissivity can 
be related to those determined for the nearby interstellar 
medium. 

Bianchi et al. (1999) have derived the FIR dust emis- 
sivity using wavelength dependences derived from FIR 
spectra of galactic emission and the SFD98 extinction map 
to normalize the tfir/Av ratio. The authors admit that 
two components of temperature could exist at low galac- 
tic latitude but they omit that dust properties could also 
vary. 

The knowledge of the extinction close to the Galactic 
plane is required to estimate precisely the properties of 
some objects, li ke the RR Lyrae ty pe stars, used as a dis- 
tance indicator, ^tutz ct al. (199E ) used the SFD98 map 
to deredden RR Lyrae type stars in Baade's window, i.e. 
3° from the Galactic center. 

5. Conclusion 

We have presented an extinction map of the Polaris high 
latitude cloud derived from star counts. This extinction 
map is compared with that produced by SFD98 based on 
the IRAS and DIRBE FIR sky maps. Within the Polaris 
cirrus cloud, the SFD98 extinction value is found to be a 
factor 2 to 3 higher than the star count values. 

We propose to relate the extinction discrepancy to 
variations in the tfir/Av of interstellar dust from the 
warm and cold emission components introduced by 
Lagache et al. (1998| ) on the basis of the observed dif- 
ference in small grain abundance and large grain temper- 
ature between the diffuse ISM and dense molecular gas. 
Within this interpretation, we find that the FIR emissivity 
of dust grains in the cold component is on average about 
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4.0 times (from 2.7 to 6.0) larger than it is in the warm 
component. There is a slight evidence that the evolution 
of the grain properties are gradual, in the sense that the 
tfir/-^v ratio increases when the cold dust temperature 
decreases. The increase in the FIR opacity could be trac- 
ing a growth of dust grains by coagulation. Our work is 
thus suggesting that the size/porosity of the large dust 
grains evolve within low opacity molecular cloud, such as 
the Polaris Flare. 

SFD98 have converted the FIR opacity in extinction 
assuming a single emission temperature for each line of 
sight and tfir/Av ratio. Our work questions the validity 
of this assumption in all regions with significant cold emis- 
sion (Fig. A few independent studies in other molecu- 
lar clouds support this statement. The SFD98 extinction 
could be overestimating the true extinction in regions with 
cold emission that are found even at high Galactic latitude 
and for extinction as low as Ay <^ 0.5. Variations in the 
dust emissivity preclude the conversion from FIR optical 
depth to Av with one single factor valid over the whole 
sky. Low latitude regions {b < 10°) and a small fraction 
of the high latitude sky require special attention. 

Acknowledgements. L. Cambresy acknowledges partial sup- 
port from the Lavoisier grant of the French Ministry of Foreign 
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